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ABSTRACT 

The matter density field exhibits a nearly lognormal probability density distribution (PDF) after 
entering into the nonlinear regime. Recently, it has been shown that the shape of the power spectrum 
of a logarithmically transformed density field is very close to the linear density power spectrum, 
motivating an analytic study of it. In this paper, we develop cosmological perturbation theory for the 
power spectrum of this field. Our formalism is developed in the context of renormalized perturbation 
theory, which helps to regulate the convergence behavior of the perturbation series, and of the Taylor- 
series expansion we use of the logarithmic mapping. This approach allows us to handle the critical 
issue of density smoothing in a straightforward way. We also compare our perturbative results with 
simulation measurements. 

Subject headings: cosmology: theory — large-scale structure of universe 
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1. INTRODUCTION 

The large-scale structure of the Universe is typically quantified with the overdensity field. This is for good reason; 
the overdensity field (5(x) = p{x.)/p — 1 (with p the density) seems to be essentially a Gaussian random field at early 
epochs, and also d sources the gravitational potential in the cosmological Poisson equation. However, when fluctuations 
are large at late times and small scales, 6 acquires some inconvenient features. It cannot go below —1, but becomes 
arbitrarily large in non-linear structures; this in itself forces non-Gaussianities when fluctuations become large. 

On small scales and late times, the probability density distribution (PDF) of the de nsity field has a ro u ghly lognormal 
form. This was noticed regarding the galaxy density field by e.g. iHubbid |l934l) and LHamilton] (|1985l ). IColes fc JonesI 
([199 LI explored the theoretical origins of the roughly lognormal nature of the matter density field. A lognormal model 
for the density field shows many features of the true nonlinear density field, such as the constraint that S cannot go 
below —1. As they showed, the continuity equation can be simplified by considering a logarithmic transformation of 
the density field A(x, r) = ln[l + (5(x, r)], 

dA{^,T) d ln[l + J(x,r)] ^ , , ^. , 

^ = ^ = -V-u(x,r) = -6'(x,t). (1) 

Here the Lagrangian total time derivative d/dr = d/dr + u-V, 6 ^ V • v, r is the conformal time, and u is the peculiar 
velocity. This equation is about as simple as one could hope for, although some complexity is hidden in the mismatch 
between the Lagrangian time derivative and the Eulerian divergence. In the limit that 9 grows according to linear 
theory, given Gaussian initial conditions, this predicts an exactly lognormal PDF for particle overden sities. Indeed, 
the v elocity field remains close to Gaussian even when the density field enters the nonlinear regime (e.g. iKofman et al.l 
Il994[ ). giving a nearly lognormal PDF for i5. 

A is also a natural density variable in Schrodinger perturbation theory (jSzapudi fc Kaiseill2003[ ). in which tree- level 
perturbation theory for A corresponds to partial loop summation to infinite order in S. This apparently captures the 
most impor tant terms, expla-ining most of the Eulerian perturbation theory loop-level variance in S. 

Recently. iNevrinck et al.l (l2009fl an a lyzed the power spectrum of A in the high-resolution Millennium simulation 
(ISpringel et al.ll2005l ). iNevrinck et al.l ([2009h found the nonlinearities in both the shape and covariance matrix of the 
A power spectrum to be dramatically reduced compared to the conventional S power spectrum. The shape of the A 
power spectrum seems intriguingly close to linear theory down to fc ^ 1 /iMpc~^, with S smoothed onto a 2-h~^ Mpc 
grid. There are a few possible reasons for this near-linearity of the A power spectrum. It could be understandable with 
perturbation theory, a possibility which we begin to explore in this paper. It could be linked with the suppression of 
density peaks in regions that have undergone shell crossing. Or, it could just be a coincidence of particular aspects of 
a AGDM-like power spectrum, e.g. the small-scale slope. 

In this paper, we make the first perturbative prediction for the A power spectrum. There are reasons to think that 
perturbation theory for A might work even better than for S. At small scales, the variance in A is much smaller than 
in S in the nonlinear regime. Also, it is tempting to try to exploit the simplicity of Eq. ([T]) in a Lagrangian approach. 

However, working with A also presents some challenges. For the S field, smoothing (onto a grid, for example) has 
fairly trivial effects; at large scales, it leaves phases intact and dampens amplitudes ne gligibly. But smooth ing is 
crucial to consider for A, giving different large-scale biases for different smoothing scales (|Neyrinck et al.l [20091 ). It is 
quite difficult to model final-conditions Eulerian smoothing using Lagrangian variables; this is why we use an Eulerian 
approach. For this, we use a Taylor series expansion of the logarithm; the slow convergence of this alternating series 
in itself portends low accuracy at highly nonlinear scales where 6 is large. 

Despite these challenges, the perturbative calculations we present below achieve good agreement with simulation 
measurements. In the next section, we review the basic framework of Eulerian dynamics and discuss the logarithmic 
transformation. We calculate the power spectrum of the A field in standard perturbation theory in section 3, and in 
renormalized perturbation theory in section 4. In section 5, we present our calculation and compare it with simulation, 
and we conclude in section 6. 

2. EULERIAN DYNAMICS AND THE LOGARITHMIC TRANSFORMATION 

The gravitational dynamics of a pressureless fluid before shell crossing is governed by the continuity, Euler, and 
Poisson equations. 

^^|l^ + V.[(l + <5(x,r)) u(x,r)]^0, 
9u(x, r) 



dr 



+ 'H(r)u(x, r) + u(x, r) • Vu(x, r) = -V<I>jv(x, r) 

vHn{^,t) = ^n„,{T)n^{T)S{^,T). (2) 

Here, T-L = d\na{T)/d\nT is the Hubble expansion rate, a(r) is the scale factor, ^m{T) is the ratio of matter density 
to critical density, and $7v(x, r) is the Newtonian potential. In principle, perturbation theory of the log-density fleld 
is possible using the log-continuity equation, Eq. ([T]), together with the last two equations in Eq. However, we 
found it more straightforward to take smoothing into account if we expand A in terms of 5, so this is the approach 
that we take. 



2.1. Equation of Motion and Perturbation Theory 

First we briefly review results for the conventional overdensity field S. Following ICrocce fc Scoccimarrol (I2006al lb[) . 
the equation of motion in Fourier space can be written in a compact form by defining the two-component variable 



(3) 



The index a G {p, v} stands for density or velocity variables, respectively. Where numbered indices are more convenient, 
we also use a G {1, 2}, i.e. = ^fp, and — 
The equation of motion then reads 



d^^aiK V) + ^ab^biK V) = 7abc(k, ki , kj) *b (ki , 7?) *c (ka , ??) , 



(4) 



with the convention that repeated Fourier arguments are integrated over. The time rj ~ In a(T) in a Einstein-de Sitter 
(EdS) universe. Here the constant matrix 



^nh — 



-1 

-3/2 1/2 



(5) 



derived for an EdS universe, is still applicable in other cosmologies with negligible corrections to the coefheients, using 
77 — \nD{T) (with D the growth factor). 
The symmetrized vertex matrix 'jabc is given by 



7222(k, ki,k2)=(5D(k- ki - ka) 



7i2i(k, ki,k2)=(5_D(k- ki - ka) 



|ki+k2|2(ki.k2) 



(ki 



ki 



2fc2 



(6) 



7ii2(k, ki, ka) = 7i2i(k, k2, ki), and 7 = otherwise. 
Then the formal integral solution to Eq. (U) can be derived as 



*a(k,?7) =.9afc(?7)0b(k) + / d??' ,ga6(f?-?7') 7bcd(k,ki,k2) *c(ki,?7')*d(k2,f7') (7) 

Jo 

where 0a (k) denotes the initial condition 0a (k) = ^^(k, = 0), and the linear propagator gabiv) is given by 





'32' 


e-3')/2 


■-2 2 " 




32 
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(8) 

In the following, we adopt growing-mode initial conditions 0a(k) = (5o(k)wa. Here Ua — (1, 1) is a unit vector. 
A perturbative solution to Eq. ([7]) is then obtained by expanding in terms of initial fields 

oo 

*a(k,77)-E*^"^(''''^) 

n=l 

*i")(k,77)= j d3qi...„ Jc(k-qi...„)J-^;\\..,Jqi,-.-,q„;77) 06, (qi) • • • 06„ (q„) (9) 

where d^qi...„ is short for d^qi • • • d'^qn, and cx^...n. denotes qi + ••• + ^ n- The kernels J"^"^ are fully symmetric 
functions of the wave vectors. As shown in ICrocce fc Scoccimarrol (|2006bD . they can be obtained in terms of gab and 
7abc recursively. 



Jl"^(ki, • • • , k„; 77)5D(k - ki...„) 



J2m=i lo ^V'dabiv - v') 7bcd(k,ki...,„,km+i...„) J'i™'(ki...m;?7') 



(10) 



symmetrized 



For n = 1, J-a^\ri) — gab{il)'^b- It should be noted that in this f ormalism, the kernel depends on the time ?7, since 
it includes subleading terms in (jCrocce fc Scoccimarrol l2006aD . If one only considers the fastest-growing mode, 

Tci^^T]) equals the well-known FT kernel ex.p{nr]){F^"\G^"^. 

2.2. Logarithmic Transformation of the Field 
In general, one can define a new vector $a as a nonlinear transformation T of *a, ^a(k, r]) = ($p(k, r]), $„(k, rj)) = 
T[\I'a(k, 77)] . In this paper, we consider a logarithmic transformation of the density field in configuration space, with 
$p(x, 77) = ln[l -|- ^'p(x, 77)], leaving the velocity field unchanged, i.e. $„(x, 77) = \l/^(x, 77). 
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$p(k, 77) = A(k, 7,) = /: [A(x, r;)] = J2 * • • • * * J (k, r;) 



ri=l 



(-1)"+! f 

= 22 / c^^qi-n '5i5(k-qi...„)*p(qi,?7)---*p(q„,77), (11) 

ri=l 

where £ represents a Fourier transform, and * denotes a convolution of overdensity fields. In the following, we will 
interchangeably use ^'p and S, $p and A when the meaning should be clear from context. Substituting the perturbative 
series expression of 4'p(q, 77) Eq. (|9]) into Eq. (|TT|) . wc get the formal expression of $p(k, -q) in terms of the initial density 
field 

°° f— 1')"+-'- f r 

$p(k, ry) = ^ / d3ki...„ J,,(k - ki...„) / ^'Pi -™i • • ■ ^'Pi-™. 

n— 1 mi , ■ ■■.ITLn 

5d{\^i - p}...™ J • ■ • (51) (k„ - Pi...™ JJ^^I^I^.^.bi^^ (p}, • • • , p^^ ; 77) • • • ^^f!!^{}.b^^ (p'l , • • • , P™„ ; ?7) 

06} (P}) • • • '^61.^ (pL, ) X • • • X (pl^) . . . ^ (p;; J. (12) 

With this definition, it is possible to express $p(k) in a diagrammatic way. As depicted in Fig.([T]), each diagram con- 
tains two levels of interactions. The thick solid n-branch tree represents the nonlinear convolution, i.e. the integration 
in Eq. ()lip . and their interaction vertex, shown as a solid square, carries a Taylor expansion coefficient (— l)""*"^//!, 
where n is the number of branches. Each branch right after the grey square represents one nonlinearly evolved density 
contrast \l/p(k„, 77) at time 77, which is followed by another m thin branches, to express its gravitational nonlinearity in 
terms of initial fields and perturbative kernels . At the end of each branch, small open circles represent the initial 
conditions </>. In short, in this diagramn iatic representation, we .j ust combine n conventional diagrams representing ^ p 
(e.g. diagrams at the right of Fig. 2 in iBernardeau et all (|2008l )) and glue them together at the solid square. When 
expanding $p(k, 77) to different orders, one should also include the information about the number of topologically 
equivalent diagrams, which essentially comes from the multinomial coefficient of (^l/p)". The final field 3>p(k, 77) is then 
a summation of all possible diagrams. 




Figure 1. Diagrammatic representation of <I>p(k, r;) up to third order. Every single branch right aft er the solid squar e repr esents one 
nonlinearly evolved density contrast ^! p and is i dent ical t o, fo r example, diagrams at the right of Fig. 2 in IBernardeau et al.l II2008I ). Explicit 
expressions for each diagram are shown in Eq. I I15II and 1 II6II . 



The power spectrum P$ (^,77) is defined as 

P$^(fc,7y)Jz5(k + k') = PA(fc,77),5c(k + k') = ($p(k,ry) $p(k',77)) (13) 

In this paper, we assume Gaussian initial conditions, with all the statistical information encoded in the initial power 
spectrum 

(0,(k)0fc(k')) ^5d{^ + ^')Pab{k), (14) 

where Pab{k) — UaUbPQ{k) for growing-mode initial conditions. Diagrammatically, the ensemble average is obtained 
by gluing two open circles together to form the symbol ®, which represents the initial power spectrum Pabik). 

Before proceeding, let's examine the definition of $p carefully. The overdensity can be idealized as a continuous 
field. In practice, however, estimating a density field from points (galaxies or simulation particles) requires a choice of 
a finite resolution, effectively a smoothing length. Furthermore, since the value of density field ^'p(x) strongly depends 
on the scale of the smoothing, as we will see, the resultant field as well as its statistical properties also vary according 
to the smoothing procedure adopted. In the following, we will consider the density field ^'p smoothed by a spherical 
top-hat window function with radius R. 



3. STANDARD PERTURBATION THEORY 

Given the perturbation expansion of $p(k, r/) as well as the dynamics of the density field, it is straightforward to 
calculate the power spectrum of a logarithmically transformed field within the standard Eulerian perturbation theory 
(SPT). Up to third order, it follows from Eq. ([TT]) and Eq. ^ that 

*p(fc) = + -fi^Hk) + <i>(3)(fc) + • • • 

(fc)V--- (15) 



1 r 
3 



P P P 



As discussed previously, a smoothing procedure is usually involved in practice. Applying a spherically symmetric 
kernel W{kR) of characteristic scale R, we get 

$W(fc) = 5'W(fc)H/(fci?) 

$(^)(fc) = vI/(2)(fc)W^(fci?)-i I d'p vI.W(p)xpW(|k-p|) W{pR)W{\k-p\R) 

^f\k)^-^f\k)W{kR)- j d^p *(,2)(p)^a)(|k_p|) wipR)Wi\k-p\R) 

+ 1 J rf'Pi2 *W(Pi)*|,'Hp2)*^,'nik - P12I) M^(pii?)W^(p2i?)W^(|k - Pi2\R) (16) 

These expressions are shown diagrammatically in Fig. [T] ($p^'' and <i>p^'' in the top row, and $p'^'' in the bottom). Note 
that the smoothing is taken over the evolved nonlinear density field 4'p. The power spectrum Pa(^) is then obtained 
by calculating the series 

PA{k) = Pi''\k) + 2 Pi")(fc) + P^f\k) + ■■■ (17) 

where P'^^\k) = ((i>p''(k)<i>p"'''(k)). It is clear that the tree-level power spectrum PA{k^rf) equals exactly the one of the 
standard overdensity field 

Pi''\k,7j) = P^^'\k,7j)W\kR) = D^(7j)W^(kR)Po{k), (18) 

where D(r]) is the linear growth function, and Poik) is the initial power spectrum. In the following, we will also denote 
the linear overdensity power spectrum Pg^^\k, rj) as Phik^ 77). To see the crucial role played by the smoothing window 
function, let us calculate 2P^^\k,ri). 

2Pi'''\k,rj)^2P^'^\k,f])W\kR) + J d^p,2 (<5o(k),5o(Pi)'5o(P2)^o(-k-pi2)) 



^ D\7j) W{kR)W{piR)Wip2R)W{\k + pulR) - 2T^^\p2, -k - P12; r?) 
xW{kR)W{piR)W{\k + pi\R) 



(19) 



where 60 is the initial fluctuation. This can be further simplified by performing the angular part of the integra- 
tion. Consid ering only the fastest-growing mode and assuming spherical top-hat smoothing (see [Appendix A| and 
[Appendix B| for more details), Eq. ([T9)) becomes 



2Pi"'(fc,,y) = 



2Pr(M)~(^ + i^)4/'.(M) 

2 



WlikR) 



21 3 dlni? 

^4 PL{k,v) kR Wl^{kR)WthikR). (20) 
where the variance of the smoothed density fluctuation cr\ — j cJ'^q PLiq^v) ^th(9-^)' ^^"^ derivative 

4^^ = ^ / '^'^ '^^ WUqR)W,^iqR). (21) 

Here we have used the geometric properties of the spherical top-hat window function derived bv lBernardeaul (|1994ai rbl. 
[1996,). 
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(13) (13) 

Since Pj^ is proportional to the linear power spectrum Pl, and it is well-known that {Pg / PL){k ^ 0) — > 0, the 
power spectrum at large scales is biased with respect to the density power spectrum 



Without the smoothing term, the bias in Eq. (|22|) would equal 1 — 26/2 1 a%,. This exp r ession accords with the 
first-order perturbative bias in Eq. (A2) (after expanding the cumulants) of INevrinck et all (j2009[ ): it is also plotted 
as a function of redshift in Fig. 4 of that paper. The bias becomes unphysically negative at z < 4 for 2-h~'^ Mpc cells. 
This smoothing term increases the predicted bias, extending its range of validity to significantly smaller smoothing 
scales. However, we found that it still goes pathologically negative at small smoothing scales. This issue is resolved in 
a renormalizcd approach, which we describe in the next section. 

4. RENORMALIZED PERTURBATION THEORY 

As the density fluctuations evolve into the non-linear regime at late times, the validity of standard perturbation theory 
breaks down. Loop contributions become ill-behaved and the convergence of perturbation series gets out of control. 
Several different approaches beyond SPT have been proposed recently, e.g. the renormalizcd perturbation theory (RPT) 
([Croccc & Scocci marro 2006a.b; Bcrn ardcau et al. 200^, the Lagrangian resummation theory ( Mat subara 2008|) , the 
closure theory (jTaruva fc Hiramatsul 120081 ). and the time r enormalization group t heory (T RG) ()Pietronil 120081) . In 
this section, after a brief review of the RPT introduced by iCrocce fc Scoccimarro! (|2006al lbf) , we will show that it is 
possible to construct the perturbation series for the transformed density field based on the non-linear propagator of 
the density. 

4.1. Renormalized Perturbation Theory of the Overdensity Field 

The crucial step of renormalized perturbation theory is to define the generalized growth factor, known as the 
propagator Gab{k) 

G,,{k, ^)Snik k') = r«,,(fc, v)SDik - k') ^ (^^^) , (23) 

which effectively describes the time evolution of individual Fourier modes when non-linear mode-coupling is included. 
The propagator measures the dependence of a non-linearly evolved Fourier mode ^a(k, rj) on its initial state 0ti(k) on 
average. Intuitively, one should expect Gat to decay to zero at small scales since non-linear mode-coupling has erased 
all the information from the initial state at that k. Indeed, with the help of the Feynman diagrams introduced in 
ICrocce fc Scoccimarrol ()2006aD . the dominant contribution can be summed up explicitly in the large-k limit, and giving 
the Gaussian decay Gab{k) = exp(— A;^(7^(a — 1)^/2). In this paper, the non-linear propagator is diagrammatically 
presented as a grey circle with an incoming branch. It is a summation of infinite number of loop contributions, as 
illustrated in Fig. [2l 




Figure 2. The nonlinear propagator r^\k,ri) has an infinite number of loop contributions. 

Given above, ICrocce fc Scoccimarrol ()2006al ) was able to rewrite the nonlinear power spectrum as a summation of 
two contributions 

Psik,f]) = G^{k,r,)Po{k) + PMcik,v). (24) 

where G is the density propagator G{k,r]) = GibUb, and PMcik,r]) is the mode-coupling term. Therefore, the non- 
linear power spectrum at any k is composed of two parts. One is proportional to the initial power spectrum at the 
same k; the other comes from mode-coupli ng of other k^ As G(k) decays at small scales, more and more power 
comes from the mode-coupling contribution. iBernardeau et all (|2008f ) showed that these complicated mode-coupling 
contributions can be expressed as a summation of multi-point propagators, defined as 

ri"U....„ (ki, • • • , K;v)SD{k - ki...„) = 1 (^^^^(k*:%''j^(k„)) ' 

which is nothing but a generalization of two-point propagator Gab- When concentrating on the density contrast 
(a = p), we will denote Eq. (^5]) as r^"^^...^^^ (or ^'il'^ bi---b„) ^'^^ simplicity. Similarly, the dominant part of multi- 
point propagators can also be summed and decay into the nonlinear regime at the same rate as two-point propagator. 



iBernardeau et all (|2008D showed that a simple approximation which generahzes the /c-dependence of two-point prop- 
agators agrees with the data with acceptable accuracy, 



r!ir^(ki,---,k„;?7) 



ri^^(|ki...n|) 

t(1, tree) I 



n^' ^^^■^■^(|ki...„i) 



r|[, ■'(ki, • • • , k„; 77), 



(n > 2). 



(26) 



With Eq. (|25p . the nonlinear power spectrum can then be expressed as 

-P* ab{k,r]) = ^'r\ I d^qi...r fo(k-qi...r) r^''^^(qi, • • • , q^; 7?)r|;J'\(qi, • • • , q^; 77) Po(gi) • • • -Po(gr). 

Pmc is the summation of all terms r > 2. For density fluctuations, it reads up to one-loop order 
P5{k,ri)= \T','\k;rj)y Poik) + 2 f d^q [rf (k - q, q; r,)! ' Po(|k - q|)Po(<z) 



(27) 



(28) 

Note that Eq. ([27|) describes nothing but an alternative way of taking ensemble averages, or diagrammatic speaking, 
gluing initial states. Instead of gluing two density fields order by order, one can first construct the objects by gluing 
initial states of individual density fields with n incoming branches, known as an (n + l)-point propagator, and the 
final non-linear power spectrum is then obtained by gluing two propagators together. With this observation, we can 
generalize this renormalization method to the perturbation theory of log-transformed fields. 

4.2. Renormalized Perturbation Theory of the Log-Density Field 
In this section, we apply the renormalized perturbation theory Eq. (^51 B7| to the log-transformed field ^p. One 
can define the {n + l)-point nonlinear propagator v'"^^ of $a. 



r 



(«) 



(ki,---,k„;?7)(5_D(k-ki...„) = — 



<5"*a(k) 



bi(ki) • • •50f,„(k„) 



(29) 



For $p, we wih denote F^'^,^...,^ or F^"^\^...,^ 



Unlike its counterpart F^'\ F^'' encodes not only the nonlinear 



information of gravitational evolution, but also of the transformation. Substituting the perturbation series of $p 
Eq. pT|) into Eq. (P^ . one notices that F^-* can be expanded in terms of the number of ^ p involved. That is, 

^ * abi---b„ 



V F^"' where 
Z^p ^ * abi- -6„' wneie 



p) 

<& abi---b, 



(ki, • • • , k„; ?7)(5_D(k - ki...„) = 



n\ p 



d^qi...p (^^(k - qi...p) 
[*p(qi,r/)---*p(qp,?7); 



bA^i)^ 



(30) 



In Fig. ([3]), we depict the expansion for the two-point propagator. The dark ellipse here represents the complicated 
ensemble average in Eq. pop . Note the Taylor coefficient of the logarithmic transformation carried by each diagram. 



r 



(1) 

A 



■X 





+ 



Figure 3. Diagrammatic demonstration of the expansion of the two-point nonlinear propag ator r^^^ (Eq. |30ll. The terms on the RHS 
correspond to contributions from n = 1 and p = 1,2 and 3. 

Taking the two-point propagator as an example, now let's write down all contributions formally. First of all, it is 
clear that the one-4'p term, i.e. the first diagram in Fig. ([3]), equals G{k,ri). For the two-VE'p term, i.e. n = l,p = 2 in 
Eq. (pO)) . we can use Eq. ([9]) to further expand perturbatively 



T^l^;^\\^-ri)^-]^x2 I d\i2 <5z,(k-qi2) 



X 2y" d^qi2 Soik - qi2) ("1 + '^^ J d^Pi-n,d^Pi-n2 ^oiqi - k - pi...„J 



X 5D(q2 -pi...„J Ji"act ^■^•c„,(k,pi, 77) J"^"J^^.rf^^(p'i,---,p;^;??) 

X (0ci(Pl)---</>c„,(Pni) </>di(p'l)---</>d„,(p^^J), 



(»2) 



(31) 
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where ni, 712 are expansion indices of two density fluctuations The coefficient (ni + f ) comes from the functional 
derivative S'i>p/6(j)a- Note that we have integrated out the variable Pm+i in the last equality together with the Dirac 
delta function SoiPni+i — k), which resulted from the functional derivative. 

By Wick's theorem, the joint ensemble average in the second equality can be written in terms of combinations of 
two-point correlations. Each of these terms can be labeled by three indices ri, r2, t, where ni ~ 2ri +t,n2 = 2r2 + t, 
i.e. we classify all pairs into three categories: ri pairs within the first field 0c;(pi), (1 < i < rii); r2 pairs within the 
second field (j)dj{p'j), (1 < j < ^12); and t pairs in between. The renormalization is then achieved by realizing that a 
pre-summation of ri and r2 gives rise to one {t + 2)-point and another {t + l)-point propagator. Therefore, 



r^',''(k;77) = -^(n + i)! J dV-n r("+^^(k, pi, • • • , p„; 7y)r(" Vpi, ' • • , -p™; ^) ^o^i) • • • Pob„), 



where we have introduced the notation 



ri"^(pi, • • • ^Pn;??) = r^"p(pi,---,Pn;f7) = r^' 



— -r(") 



(") 



pci ■ 



(pi,---,p„;?7) 



(32) 



(33) 



and used the definition of T^^ in terms of the perturbation kernel. The coefficient n\ then comes from all possible 
ways of matching n initial states with another group of n initial states. It is not hard to notice the similarity between 
Eq. ([5^ and Eq. (I?7|) . since the derivations are nearly identical. Yet, there are also several differences. For the RPT 
description of the matter power spectrum, both propagators share the same order; otherwise, it would be impossible 
to match the pair. For the same reason, the order of one propagator in Eq. (j32]) is always greater than the other by 
one, since we have to select one branch out before taking the average. To conclude, we are able to construct from 
terms involving F^. 

Similarly, for the three-^p term, we have 

F^^'/)(k;,7)= J2 fn^n^.n,,/ ^'pi- ■ dVi- . dVi' . .„3 r("^+"^+^) (k; pi , • • • , p„, ; p'^ , • • • , p'„^ ; 7?) 



X F 



("2 +"3) 



(-P'l,---,- 



("3+"!) 



(-P'l'r 



Pi, 



X PoiPl) ■ ■ ■ PoiPn,) Poip'l) ■ ■ ■ Poip'nJ Mp'I) ' ' ' Po^Pn,)- 

where, for each contribution 



(34) 



: (ni +712 + 1) 



ni + 71,2 



"2 + 
"-2 



"-3 + "1 
"3 



ni!ri2!'^3' 



(35) 



In this case, in order to describe all configurations, we need 3 indices to label the internal pairing and another 3 indices 



(71,1,71,2,71.3) for pairing in between. The coefficient of each contribution gn^^ 



can be understood as the number 



of equivalent ways of gluing pairs given the triplet ric = (711,7x2,713). The number of initial states for each 5'p is 
Hi, = (77,1 +772 + 1, 772 + ?T-3, + 77i). After selecting one initial state from the first \E'p (with 711 + 772 + 1 possibilities), 
one has to count the number of choosing ric out of n„, which equals ("^^j"^) ("^^"^) ("^^"^) • The formula can be 
further generalized to arbitrary order 77, 



.(1, n) 
A a 



(k;r,), 



n '^P^-u 

l<'i<j<n 



n [p^(p''^] 

l<'i<j<n 



tij 



-p(ti+i;/i . 11 



Pi 



nn 



(36) 



where the index tij denotes the number of connections between the i-th and the j-th density field ^'p (e.g. 
{^12,^13,^14,^23,^24,^34} for 77 = 4), so the total number of the indices is then 77(77 — l)/2. The coefficient 
then 

9{u,} 



,(1, ") 



(-1)"+^ X {ti + 1) 



n 



til ' ' ' ti 



n (*»^)' 

l<i<j<n 



(37) 



Here, ti = tim is the total number of connections linked between the 7-th density field 5'p and others, with tij = tji. 
Eq. p7p expresses the products of multinomial coefficients of choosing {tn, • • • , from ti for each density field times, 
all possible permutations within tij for pair matching. Each connection carries a momentum p^"*, (1 < I < tij), which 
characterize the Z-th connection between j-th ^p and j-th ^fp. The integration is taken over momenta of all possible 

connections among 77 different density fields. Because of the ensemble average Eq. (ITi|) . we have p]-'.^ = — Pf*;. And 
we also used the shorthand notation [Po{p^-')Y'^ for Po{pY) ■ ■ ■ PoiPu )■ 
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2x < c f^ < J 3 — <- 

k — q 




Figure 4. Nonlinear propagator r^'(fc) of ln(l + S) up to one-loop order, corresponding to contributions in Eq. II38I I and I I39I I. 



Figure 5. Each diagram contains an infinite number of loop contributions. 



Now we can expand in terms of the number of initial power spectra entering the calculation. For the tree level, 
it simply reads as 



r(^';'™)(k;,y)=r«(k;ry) W{kK) 



(38) 



At one-loop, from Eq. (j36p . the contribution would be nonzero only when n < 3, since otherwise there would exist at 

least one Ts with zeroth order, which would vanish because F^"'' = {^p) = 0. For n — 2, ti2 — 1, the coefficient equals 
—2. For n = 3, i23 = 1, ^12 = ^13 = 0, the coefhcient equals 1. So we have 



t(1, 1 — loop)^ 
A a 



(k;r;)= / rf^qPo(9) 



WiqR) FW(q; rj)] T^ll{k; ri)W{kR) - 2 W^(gi?)M^(|k - q|i?) 



(39) 

Given Eq. p6p . we can draw every contribution diagrammatically. Starting from the diagram representing ^p"-* in 
Fig.([T|), we change all kernels into n-point propagators. After selecting one particular branch out, we glue the rest of 
the initial states (open circles) together. Every resultant topologically inequivalent diagram represents one or several 
terms in Eq. (j36p . Since all ordinary kernels have already been substituted by propagators, the ensemble average 

(gluing) is only performed among different density fields ^'p. In Fig.Q, we show all diagrams of fJ^-* up to one-loop 

order. They correspond one-to-one to Eq. (|38l) and Eq. (|39l) . We also present all two- loop diagrams of F^' in Fig. (fTO| . 
It should be emphasized that, by substituting ordinary kernels into propagators, each diagram in fact contains an 
infinite number of loop contributions at every substituting position, as shown in Fig.(IS|). 
For each diagram, there exists a straight path through the diagram carrying the same momentum k at both the 

start and end. Along this path, there is one convolution vertex, as well as one (n + l)-point propagator F^"'. For 

n = 1 (e.g. the third diagram in Fig.Q), this contribution recovers the same fc-dependence as F^^''(fe), rescaled by a 
constant from loop integration. When n > 1, Eq. (|26p suggests a similar damping of k given pi • • • Pn-i- Meanwhile 
every propagator associates with a smoothing window function W{kR). Therefore one should expect that, in the 

large-fc limit, F^^ will decay as a combined effect of a Gaussian damping of F^^^ and the smoothing window function. 
Furthermore, as we will see in the next section, the presence of the propagator within the loop helps to regulate the 
convergence of the numerical calculation. 
One can also derive the three-point propagator similarly, except that two distinct contributions have to be taken 



into account because of the second derivative in the definition of F 



(2) 



F^\;)(ki,k2;^): 



2n 



X n 



+ (n- 1) 



qi - n ^D(k - qi...„ 
(5*p(qi,77) (5*p(q2,?7) 



(52*p(qi,77) 



^'p(q2,?7)---*p(q„,77) 



*p(q3,?7)---*p(qn,?/) 



(40) 



As shown in Eq. (|40p . the first term takes the joint average of the products of a second derivative with n — 1 density 
fields. Diagrammatically speaking, this category includes the first, third and sixth diagram in Fig.Q, i.e. both 
incoming branches come from one single density field ^'p. The rest of the diagrams then correspond to the second 
term of Eq. (|40p. where two branches originate from two different density fields. Both terms can be written in the 
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Figure 6. Three-point nonlinear propagator F^' up to one-loop order, corresponding to contributions in Eq. I|43|l and l|44|l . 

same form of Eq. (j36p , 

(2, n) ^(ti+2). 



{tij} l<i<j<n 



X ^^)(p-,„^,...,pr..„J + ffjy r(*^^)(ki;p}i,^^,...,p-,^j 

X ri*f ^)(k2; p?^.,,, , • • • , ■ ■ • ri*")(p^.i,„^, . . . , pr,„ j| n Wt' 

l<i<j<n 

(41) 

The difference between the two terms can be clearly seen from their orders of their propagators. With the same 



labeling system, the first gives 
g coefficients equal 



5 •••-^5 



" , while the second gives T 



(tl + l)T.(t2+l) T.(t„) 



rV" . Meanwhile, the two 



n 

l<i<n 



' ' ■ ti' 



n (*»^)' 



l<2<_7<n 



/I \n+l 
J2, n) 



n 

l<i<n 



n 



l<z<j<n 



(42) 



The first contribution g)^^.^' gives n(ti -I- 2){ti + 1), while the second term gives n{n — l){ti + l){t2 + 1). 



(2) 

At tree level, there are two diagrams: one is , from the single- vE'p contribution; and the other is the two-\E'p term 
with ti2 = for the second term of Eq. (HT|) . 

r^^'„r)(p, k - p; ,y) = rflip, k-p;rj) W{kR) - (p; ^)rW (k - p; rj) WipR)Wi\k - p\R). (43) 

At one-loop level, five diagrams are nonzero. 



t(2, 1-loop) 
A ab 



(p,k-p;r/)= J d^pPoiq) 



.1=1 



(2, 1-loop) 



(44) 



where 



(2, 1-loop) 



- -3 W{qR)W{\k + q\R) (-q; vK"Jk - p, p, q; 77) 



.(3) 



(2, 1-loop) 
2 " 
(2, 1-loop) 



s aty 

(2) 



-2 W^dp + q|i?)W^(|k - p - q|i?) r^-^;,(p, q; r,)rf,(k - p, -q; rj) 



4WiqR)Wi\p + q|i?)W^(|k - p\R) r^^>(q; ??)r^^i(k - p; rj)r';i{p, q; 77) 



.(2), 
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^(2, 1-loop) 



WiqR)T^P (q; ry) WikR)rYUp, k - p; ,7) 



(2) 



JC 



(2, 1-loop) 



W'iqR)W{\k-q\R)WikR) Tf\ci;v) T^^lk - q; r?)r^^;,(p; 77). 



.(1) , 



(45) 



In Eq. (|45l) . /Ci and /C2 are two-^'p contributions from the first and second terms of Eq. (|40l) respectively, IC3 and 
/C4 are three- contributions, and JC5 is the four-^'p contribution. 
At last, the non-linear power spectrum of the log-transformed field $p can be expressed as 



PA{k;v) 



r>l 



d-^qi...r SoCk- qi...r) T^^\qi, ■ ■ ■ ,qr;ri) Po{qi) ■ ■ ■ Poiqr 



5. NUMERICAL RESULTS 



(46) 



In this section, we show our numerical results. To simplify the ca l culati on, we adopt the approximation Eq. (I26p . 
which has been shown to be accurate enough in iBernardeau et al.l ()2008() . For the tree-level propagator, we only 

include the fastest-growing mode, i.e. the standard perturbation kernel F^^K Since the non-linear quantity F^"^ is 
involved throughout our formulae, all the calculation is done num erically. All results assume the fiducial concordanc e 
ACDM cosmology of simulation of the Coyote Universe suite (jHeitmann et al.l 120101 120091 : iLawrence et al.l I2010D , 
with {V.mh'^, ^bh^, w, as, h) — (0.1296, 0.0224, 0.97, -1, 0.8, 0.72). We use these parameters because we compare 
our perturbative results to m easurements froni this simulation. The linear power spectrum is calculated with the 
public Boltzmann code CAMB ([Lewis et al.ll200(il) , and the numerical integration is performed with the CERNLIB multi- 
dimensional integration routine. 




k (h/Mpc) 



Figure 7. Non-linear propagator over linear growth rate at a = 1 and a = 0.5. The dashed line shows the two-point propagator of the 
overdensity field T^^'. Solid lines show the propagator for the log-transformed field f'^', with top-hat smoothing radius R = 20 h^^ Mpc, 
15 Mpc, 10 Mpc and 5h~^ Mpc, from top to bottom. 



In Fig. ([7]), we show the two-point propagator fJ^ (fc) up to two- loop order. As expected from the analysis of SPT, the 
propagator of the A field approaches a constant less than one at large scales. The value of this large-scale bias, which 
also encodes the statistical information of the density fluctuation fEa l^^ . strongly depends on the smoothing process 
adopted before the transformation. For the smoothing radius i? — > 00, effectively no loop integration contributes 

to F^^(fc), therefore the bias approaches to 1 from the tree-level result, Eq. ([38]). When the smoothing radius goes 

smaller, since F^"'(fc) that appears in the loop integration also decays at large k, the effect of the smoothing becomes 
less and less important until the bias freezes at some value, where higher-loop contributions are significant. In this 

sense, the presence of F^"-* (k) helps to regulate the convergence of the perturbative series. 

At small scales, T^^\k) is dampened both by nonlinearities and by smoothing, while F^^''(fc) has no additional 

smoothing imposed. This is why fJ^-* (fc) decays at larger scales than F^^-* (k) for large smoothing scales. The top solid 

line shows F^^(fc) with a smoothing scale R — 20 /i^^ Mpc, which reduces to 0.5 around k ^ 0.1/iMpc^^. As the 
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smoothing radius becomes smaller, the curve approaches Tg (k) at large k. For R — 5h ^ Mpc (the bottom solid 

line), the damping of r^''(fc) is similar to that of Tg^\k). 

If the large-scale bias is divided out, boosting the A propagators to line up for small fc, the A propagator may slightly 
exceed the S propagator on small scales, if the smoothing scale is sufficiently small. We have found this to be the case 
in preliminary simulation measurements. Disappointingly, this implies that the logarithmic transform by itself does 
not help appreciably to reconstruct mode-by-mode initial phases and amplitudes. But the similarity in damping is not 
surprising, given that e.g. bulk displacements from the initial conditions affect both fields. 




1 1 1 1 1 1 — r- 
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1 1 1 1 1 1 1 1 1 1 1 

□□^ □ p rwi n n n m , n-r. , " 




, 1 



0.01 



0.1 



k (h/Mpc) 

Figure 8. The nonlinear power spectrum of A - 



■ ln(l -I- <5) at a = 1 and a - 



k (h/Mpc) 

■ 0.7, with top-hat smoothing radius R - 



3.7 h ^ Mpc. Upper 



panels: Triangles and squares are the <5 and A power spectra measured from a simulation. Solid lines include both the F*^' and F*^' terms, 

and dotted lines include only the f'J' term. Dashed lines are the linear power spectrum, damped at small scales with the NGP window 
function, to show the degree of nonunearity of the S power spectrum. Note that we did not compute a perturbative prediction for the 5 
power spectrum in this paper. Lower panels: the nonlinear power spectrum of A divided by the linear density power spectrum. Dotted 
horizontal lines indicate the value of the large scale bias. In these figures, we use a different Fourier-transform convention; compared to 
previous sections, here power spectra are multiplied by (27r)^. 



In Fig.®, we illustrate the non-linear power spectra of both PA{k) and Ps{k) at a = 0.7 and a = 1. The squares 
and triangles are the S and A power spectra measured from Coyote Universe simulation 0, which has 1024^ particles in 
a cubic 1300 Mpc 1 h^^ Gpc) box. Its resolution is high enough that the power spectrum is accurate to sub-percent 
level down to l/iMpc~^. We measure both power spectra on 128'^-cell grids , using nearest-grid-p oint (NGP) density 
assignment. Even for the A field, the shot noise is negligible at this cell size (jNevrinck et al.ll2011[) . We do not correct 
either measured power spectrum for the NGP pixel window function. In fact, such a correction would be ill-defined 
for the A field. The dashed line shows the linear power spectrum, after applying the NGP window function. Note 
that we did not compute a perturbative prediction of in this paper that would allow a fair shape comparison to our 
perturbative prediction of Pa- 

In Fig.®, we show the same results at rcdshift a = 0.5. The calculation is quite accurate before it seems to break 
down a t fc ^ 0.31/iMpc~^. In the inside panel, we show the BAO wiggles of both PA{k) and linear Ps{k), dividing 
out the lEisenstein k Hul ()1998D ' no-wiggle' power spectrum. 

The PT result of PA{k) (solid line) includes contributions from F^^ as well as F^-*. Since the numerical integration 

involving F^'' becomes extremely time-consuming, we only calculate it up to one-loop order, while F^^ is done to 
two-loop order. The agreement with the data measured from simulation is good, given the several approximations we 
have made. For a = 1, the calculation breaks down after k > 0.2 hMpc^^ . Also, the discrepancy with data between 

k = 0.1 and k = 0.2/iMpc~^ peaks at over 10%. This is where F^"* becomes dominant, but is only calculated up to 

one-loop order. The dotted line presents the F^^ contribution. For higher redshifts, the differences between data and 
theory become smaller. At a = 0.7, the prediction breaks down around k ^ 0.28/iMpc~^. In the lower panel of each 
figure, we also plot the ratio between PA{k) and the linear power spectrum, after applying the pixel window function. 



6. CONCLUSION 
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k (h/Mpc) 

Figure 9. The nonlinear power spectrum of A = ln(l + 5), at a = 0.5, with top-hat smoothing radius R = 3.7ft.~^Mpc. Li nestyles and 
symbols the same as in Fig.|8] Inset: The solid line shows the BAO wiggles of P^(fc), dividing out the Eisens tein &: HuT lllQQg'l 'no-wiggle' 
power spectrum. The dot-dashed line presents the wigglincss of the linear Pg{k), multiplied by a constant to line up with PA{k) at large 
scales. 



In this paper, we developed the cosmological perturbation theory for the log-density field A = ln(l + S). In the 
context of standard perturbation theory, we showed how different smoothing scales can alter the amplitude of the 
A power spectrum at large scales. With the help of multi-point propagators developed in renormalized perturbation 
theory, we constructed the building blocks for the power spectrum of the A field. In our formalism, each diagram 
effectively includes infinitely many loop contributions. The Gaussian damping of the propagator for the A field ensures 
that the convergence of the series is well-controlled. 

We found that except for the large-scale bias, this A propagator, which quantifies the memory of initial conditions 
as a function of scale, is similar to the 6 propagator, damping at a similar length scale. This means that the memory 
of mode-by-mode initial phases and amplitudes in the A field is not much better than in the 6 field. 

Even with several approximations, our PT calculation for Pa achieves good agreement with simulation measurements. 
However, further work is necessary to obtain results with sufficient precision to resolve the question of whether the 
much-reduced nonlinearity in the shape of Pa compared to Ps is understandable perturbatively. 

We thank Katrin Hcitmann and Adrian Pope for help accessing the Coyote Universe simulations. XW, MN and AS 
are grateful for support from the W.M. Keck and the Gordon and Betty Moore Foundations, and IS from NASA grants 
NNG06GE71G and NNX10AD53G, and from the Polanyi Program of the Hungarian National Office for Research and 
Technology (NKTH). 

APPENDIX 
Appendix A. ONE-LOOP SPT 
In this appendix, we present the one- loop result of PA{k) in the context of standard perturbation theory. Starting 
from Eq. ([T9l) . if one only considers the fastest-growing mode, 2P^^'^^ becomes 



4 j d^p p^(p,,y)i^(2)(k,p)VF(pi?)VF(|k + p|i?) 



PL{k,Ti)W{kR), (Al) 



where we have already expressed the four-point joint average in terms of the power spectrum. i^(^)(ki,k2) = ^ + 
5// -f +1 is the second order perturbation kernel, where /x = It can be further simplified by integrating 
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out the angular part of the third term of Eq. (jAip 



W^(|k + p|i?), 



(A2) 



where fi is the cosine of the angle between the vectors k and p. To proceed, let us consider the spherical top-hat 
window function, which in Fourier space is defined as 



WthikR) 



{kR)3 



[sin(fci?) - kRcos{kR)] . 



(A3) 



For a given vector k, one can prove the following integration identities (|BernardeaTjll994ai rbl. 11996( 1 : 

J ^(1 - M')M^t/.(|k + p\R) = ^ Wth{kR)Wth{pR); (A4) 

Wth{kR)Wth{pR) 



I fk p 

- + r 



W^t„(lk + p|i?) = - 



sin(ni?) , , sinffci?) 
[WthikR) pR W;^{pR) + WthipR) kR W^^{kR)] 



(A5) 



Substituting Eqs.(|A4lfclA5|) into Eq. (|A2|) . one obtains 

26 
21 



2Pi^3)(fc,r;)=2Pi")(fc,ry)1^2(fci?) - ^4 PL{k,r^)W\kR) 



j dp Anp^R Pl{p, -q) 



xWth{pR)W[^(pR) 



PL{k, v)Wl{kR) - -al [kRW^„{kR)WtH{kR)] P^k, rj), 



which reduces to Eq. (|20|) if we define the derivative of cr|. with respect to R, Eq pTjl . As for P^'^\k), one can derive 



(A6) 



P^f\k,rj)^P^^'\k,rj)W'{kR) + J d'pPL{p,v)PL{\k~p\,v) 
-F(2)(p, k - p)W{pR)W{\k - p\R)W{kR) 



-W\pR)W^i\k~p\R) 



(A7) 



Appendix B. GEOMETRIC PROPERTIES OF THE TOP-HAT WINDOW FUNCTION 

Integra tion identities Eq. (|A4p and Eq. (|A5|) are nearly the same as th e one derived in the Appendix of iBernardeaul 
(|1994al ibr). except only one angular integration is involved. Following IBernardeaul (|1994al lbl). the top- hat window 
function can be decomposed as a sum of Bessel functions 



Wth{k) = 3y/^k-^^^J3^2{k) 

°° 3 d 

Wthi\k + p|) = 3^ ^ (- + m) (M~^/'^3/2+m(fc) J3/2+m(p) J^'P 1^) , 



m=0 



(Bl) 



where Vm is the Legendre polynomial of order m. Therefore Eq. (jA4l) can be expressed as 

I ^i^i_^2^WtH{\k + p\) = 1^ f;(^+m)(Ap)-3/V3/2+™(A:)J3/2+™b) /' d^Ji[l-^?) ^V„.^,i-^,). 

'' ^ m=0 •'-^ ^ 



Since j\ dfi (1 - ^i'^) -^Vrn+ii-y) = ff^om, Eq. is proved. 

/ ^ (1 - t^^Wth{\k + p|) = 371 (fcp)-3/2 J3/2+™(fc) J3/2+™(p) = \ Wth{k)Wth{p). 

For Eq. (|A5p . by substituting x — |k + p|i?, one gets 



(B2) 



(B3) 



3 /■I'^+Pl^ X dx 
'2J\^-P\R kpR^ 



(fc2 +p2)(a.2/^2 



.2M 



4fc2p2 



sin(a;) — x cos(a;) 



WikRf^^ + WipRf''^'''^^ 



pR 



kR 



(B4) 




In this paper, we have calculated the two-point propagator up to two-loop order. Seven non-vanishing contri- 
butions are depicted in Fig. (|10l) . From these diagrams, one can write down all terms explicitly 



2-loop), 



(k;?7) - / d3pi2 Pa{pi)Po{p2) 



where 



/CI 



(1, 2-loop) 



(1, 2-loop) 
2 

(1, 2-loop) 



1^(1. 2 — loop) 



IC, 



(1, 2-loop) 



(1, 2-loop) 
6 

(1, 2-loop) 



1^(1. 2-loop) 

.i=l 



: -6 l^dpi + P2|i?)l^(|k + pi + p2\R) Ff (-P1, -p2; v)^fi{k, pi, p2; ry), 
:6 W{piR)W{p2R)W{\k + pi + p2|i?) T's'\-pi;v)r'i\-p2;v)^fi(K Pi, P2; 
:8 W{piR)Wi\pi + P2|i?)M^(|k - p2|i?) Ff (pi, P2; t;)f(3)(-pi; ^)Ff Jk, ~p2;v) 



2 W'{\pi + p2\R)w{kR) rf>{pi,p2;v)r';>{-pi,-P2;v)^'^'a{'i^;v) 



(2)/ 



6W{p,R)W'{p2R)W{\k-p,\R)T^^\p,;T^) T's'>{p2;v) rf>^{K-pi;r,) 



.(2) , 



6 WikR)WipiR)Wip2R)Wi\pi + P2\R) F^ (pi, P2; ?7)r^''(-Pi; r/)Fy^(-p2; r,)Fyi(k; rj) 



(1)/ 



.(1), 



(CI) 



3W'{piR)W'ip2R)W{kR) rf\p,;7j) Ty>{p2;v) r^'i(k;r/) 



-(1) , 



(C2) 
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